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The framework of non-extensive statistical mechanics, proposed by Tsallis, has been used to 
describe a variety of systems. The non-extensive statistical mechanics is usually introduced in a 
formal way, using the maximization of entropy. In this article we investigate the canonical ensemble 
in the non-extensive statistical mechanics using a more traditional way, by considering a small 
system interacting with a large reservoir via short-range forces. The reservoir is characterized by 
generalized entropy instead of the Boltzmann-Gibbs entropy. Assuming equal probabilities for all 
available microstates we derive the equations of the non-extensive statistical mechanics. Such a 
procedure can provide deeper insight into applicability of the non-extensive statistics. 


I. INTRODUCTION 


Complexity in natural or artificial systems may be caused by long-range interactions, long-range memory, non- 
ergodicity or multifractality. Such systems have exotic thermodynamical properties and are unusual from the point of 
view of traditional Boltzmann-Gibbs statistical mechanics. Statistical description of complex systems can be provided 
using the non-extensive statistical mechanics that generalizes the Boltzmann-Gibbs statistics [iRsj. The non-extensive 
statistical mechanics has been used to describe phenomena in various in high-energy physics Q , spin-glasses @ , cold 
atoms in optical lattices Q, trapped ions 0, anomalous diffusion dusty plasmas |T^ . low-dimensional dissipatiye 
and conservative maps in the dynamical systems [iiHil, turbulent flows [l3 |. Langevin dynamics with fluctuating 
temperature [Tsl [l^ . Concepts related to the non-extensive statistical mechanics have found applications not only in 
physics but in chemistry, biology, mathematics, economics, and informatics as well [l7l - [l^ . 

The basis of he non-extensive statistical mechanics is the generalized entropy [l[ 


Sq — feg 


9-1 


( 1 ) 


where p(/r) is the probability of finding the system in the state characterized by the parameters /r; the parameter q 
describes the non-extensiveness of the system. More generalized entropies and distribution functions are introduced 
in Refs. [13, [HI. The generalized entropy m can be written in a form similar to the Bolzmann-Gibbs entropy 


as an average of g-logarithm [l| : 


where the (/-logarithm is defined as 


Sbg = -kBy^^pju) Inp(f^) 


S<i = kBY.piT)\nq — 


hig X = 


- 1 
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( 2 ) 


(3) 


(4) 


In the limit 9^1 the (/-logarithm becomes an ordinary logarithm, thus the Boltzmann-Gibbs entropy can be obtained 
from Eq. ([T]) in the limit (/ —>■ 1 [Ei- The inverse function of the (/-logarithm is the (/-exponential function 


exp,(a;) = [1 -f (1 - q)x]l ’ , (5) 

with [a;]+ = a; if x > 0, and [x]+ = 0 otherwise. The (/-exponential and (/-logarithm appear in many equations of 
non-extensive statistical mechanics [l[- Some properties of (/-exponential and (/-logarithm are presented in Appendix I bI 
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The equilibrium of an isolated system consisting of N particles is described by the microcanonical ensemble. In the 
statistical physics it is assumed that the equilibrium in the microcanonical ensemble corresponds to equally probable 
microstates [^. therefore in the microcanonical ensemble p(/r) = l/bb, where W is the number of microstates. 
Non-extensive statistical mechanics can describe non-ergodic systems where not all available microstates can be 
reached. In this case W is the effective number of microstates, that is the number of microstates whose probability is 
not zero. When probabilities are equal, Eq. m for the generalized entropy takes the simpler form 

Sq = fee ln,j W . (6) 

In the systems with long-range interactions and long-range correlations the effective number of microstates W can 
grow not exponentially with the number of particles in the system N but slower, as a power-law of N. For such a 
systems the standard Boltzmann-Gibbs entropy ([2]) is not proportional to the number of particles in the system and 
thus is not extensive. The extensive quantity is the generalized entropy ([T]) for some value oi q ^ 1. In general, if the 
entropy Sq is proportional to the number of particles N then the number of microstates W grows as IT ~ exp^ N. 
There are two different cases: (i) g < 1 and W ~ xhe number of microstates grows as a power-law. (ii) 

q > 1 and W behaves as (1 — (g — . In this case there is a maximum value of the number of particles 

fVcrit where the number of microstates becomes infinite and thus the macroscopic limit TV —>■ oo cannot be taken. 
Because the of this complication occurring when g > 1 in this paper we consider only the case of g < 1; the value of 
g in all the equations below should be assumed to be less than 1. The case of g > 1 warrants a separate investigation 
and is outside of the scope of the present paper. 

The canonical ensemble in the non-extensive statistical mechanics is usually introduced in a formal way, starting 
from the maximization of the generalized entropy [l| . The physical assumptions appear in the maximization procedure 
in the form of constraints. However, the g-averages used for constraints are unusual from the point of view of ordinary, 
Boltzmann-Gibbs statistics. The physical justification of g-averages and escort distributions is not completely clear. 
Thus a more physically transparent method would be useful for understanding the non-extensive statistical mechanics. 
The goal of this paper is to investigate the canonical ensemble in the non-extensive statistical mechanics using a more 
traditional way, by considering a small system interacting with a large reservoir via short-range forces. Consistent 
investigation of such a situation has not been performed yet. We assume that the generalized entropy o for some 
value of g < 1 instead of the Boltzmann-Gibbs entropy is extensive for the reservoir. In addition, as in the standard 
statistical mechanics we assume equal probabilities for all available microstates of the combined system consisting of 
the small system and the reservoir. By doing so we can avoid the critique of the generalized entropy presented in 

Refs. [13, il. 

The paper is organized as follows: To make the comparison of the non-extensive statistical mechanics with the 
standard Boltzmann-Gibbs statistical mechanics easier, in Section [H] we briefly present the usual construction of the 
canonical ensemble in the standard statistical mechanics. In Section m we consider the canonical ensemble in the 
non-extensive statistical mechanics and in Section m we explore the resulting Legendre transformation structure. 
Section IVl summarizes our findings. 


II. CANONICAL ENSEMBLE IN BOLTZMANN-GIBBS STATISTICAL MECHANICS 

To highlight the differences from the non-extensive statistical mechanics, let us at first briefly review the canonical 
ensemble in the extensive Boltzmann-Gibbs statistical mechanics. The standard approach [HI . [23l| is to consider a 
composite system consisting of a system under investigation S interacting with a large reservoir R. The system S has 
energy E, the energy of the reservoir R is Er and the energy of the composite system is Etot- Due to the interaction 
the system S and the reservoir R can exchange energy. The interaction is assumed to be short-range, therefore in 
the macroscopic limit the energy of the interaction is negligible and the total energy of the composite system is 
Etot = E + Eh- For simplicity we assume that there is no exchange of the particles between the system S and the 
reservoir R. 

The number of microstates in the system S having the energy E is W{E) and the number of microstates in the 
reservoir is Here it is assumed that the numbers of microstates depend only on the energy. The short range 

interactions of the system S with the reservoir R do not significantly change the numbers of microstates and thus the 
total number of microstates in the combined system when the system S has energy E is W{E)WYi{Etot — E). The 
full number of microstates Wtot(Ftot) of the combined system is obtained summing over all available energies of the 
system S: 


W'tot(Ftot) = ^ W{E)WK{Etot - E). 
E 


( 7 ) 


3 


Introducing the entropy of the system S{E) = fee \nW{E) and the entropy of the reservoir Sb.{Eb.) = fcs In VhR(i?R) 
we can write 

WtotiEtot) = . ( 8 ) 

E 


In the extensive Boltzmann-Gibbs statistics the entropy of the reservoir 5 'r is proportional to the number of particles 
Nn in the reservoir and is macroscopically large. The sum of large exponentials can be approximated by the largest 
term, as is described in the Appendix In the statistical mechanics it is postulated that in the equilibrium the 
probability of each microstate is the same and equal 1/Wtot- Thus the most probable state of the composite system 
corresponds to the largest term in the sum ([5]). The most probable energy U of the system S corresponding to this 
largest term can be found from the condition 


This condition allows to introduce the temperature T characterizing the equilibrium: 


1 

T 


dU 


S{U) 


T^SniEtot - U). 
dEtot 


( 10 ) 


We can also consider the situation when the Boltzmann-Gibbs entropy of the system S is not necessarily extensive 
and proportional to the number of particles N in the system. If we introduce the generalized entropy as Sq{E) = 
/cb lug W{E) then the sum (I5|) becomes 

W'tot(Atot) = (11) 

E 


and the largest term is determined from the condition 

As,((/) - - C/) = 0. (12) 

Here we have used Eq. (lB4l) . We can conclude, that the temperature T is related to the generalized entropy Sq via 
the equation 


1 ijs,{u) 

T 1 + l^Sq{U) ■ 

Introducing the auxiliary g-temperature by the equation 

we get the relation 


(13) 


(14) 


(15) 


The auxiliary temperature Tq in the formulation of the non-extensive statistical mechanics based on maximization 
of entropy can appear as the inverse of the Lagrange multiplier associated with the energy constraint. Althoug Tq 
is not the physical temperature, it can have another physical meanin g. F or example, such effective temperature is 
directly related to the density of vortices in type II superconductors [^. The relation (jTSl) between the physical 
temperature T and the auxiliary temperature Tq has been proposed by various authors in Refs. Definitions of 

the temperature associated with different formulations of the non-extensive statistical mechanics have been analyzed 
in Ref. |36l |. The general requirement that composition rules of entropy and energy should satisfy to be compatible 
with zeroth law of thermodynamics has been investigated in Ref. [s^- It has been shown that formal logarithms of 
the original quantities should be additive. 

Note, that the statistics of the system S is determined by the reservoir, as we see from Eq. dUD. Therefore, it is 
more convenient to describe even such a system using the Boltzmann-Gibbs entropy. The same conclusion has been 
made in Ref. [3^ : it has been shown that physical temperature and pressure within the formalism for non-extensive 
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thermostatistics leads to expressions which coincide with those obtained by using the standard formalism of statistical 
mechanics. 

According to the central postulate of the statistical mechanics, the probability of the microstate /i ® /tr where the 
system S is in the microstate ^ and the reservoir is in the microstate /iRis 


The probability of the microstate fjL of the system S then is 

p(m) = ■ 

MR 


(16) 


(17) 


If the energy of the microstate fi is then the energy of the reservoir is Etot — E^ and the number of acceptable 
microstates of the reservoir is WR(£'tot — E^). We obtain that the probability of the microstate is equal to 


P{p) = 


WniEtot - E^) 
W{Etot) 


(18) 


Approximating the number of microstates as 


WR(£^tot - E) = 


(19) 


we obtain that the probability of the microstate of the system S is proportional to the Boltzmann factor 


P{E) = exp 



( 20 ) 


Here we used the definition of temperature (ITOl) : g^^S'R(£’tot) ~ y ■ Note, that for the justification of the exponential 
form of Eq. (12011 it is essential that the Bolzmann-Gibbs entropy of the large reservoir had very small second derivative, 

o 2 

■S'r {Etot ) ~ 0. This requirement means that the heat capacity of the reservoir 


Ck = - 


1 

T2 J^5R(£;tot) ■ 


( 21 ) 


should be very large, that is the reservoir should be a thermostat. Other possible forms of entropy (for example, the 

/TD \ 

generalized entropy Sq — fce lug ITr with q ^ 1) do not have small second derivative and thus do not lead to a good 
approximation for the probability of microstate. 

From the Boltzmann factor (PIB follows that the normalized probability of the microstate can be written as 

pip) = , ( 22 ) 

where 

Z = (23) 

M 


is the partition function. The distribution of the energy of the system E is obtained multiplying the probability 

1 o/ T^\ 

pin) by the number of microstates having energy = E. This number is equal to W{E) = ^ , therefore the 

distribution of the energy is given by 

p[E) = , (24) 

Zj 

The probability p{E) should be normalized, thus the partition function can be also written as 

E 


( 25 ) 
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In the macroscopic limit the sum of large exponentials can be approximated by the largest term, therefore 

\nZ^^SiU)--^U, (26) 

/CB KbI 

where the energy U corresponding to the largest term is obtained from the equation 

(27) 

The average energy of the system 

U = Y,EM (28) 

can be determined form the partition function: 

U = kBT^-^\nZ. (29) 

The free energy F is introduced according to the equation 

F = -kBTlnZ . (30) 

The equality 

F = U-TS (31) 

defines the average entropy S. Combining Eqs. (ESI)-® we get 

M - 1 

W~T' 

Due to approximation (1261) in the macroscopic limit the average energy U coincides with the most probable energy U 
and the average entropy S coincides with S{U). 


III. CANONICAL ENSEMBLE IN NON-EXTENSIVE STATISTICAL MECHANICS 

Now let us consider the composite system where the large reservoir R is not described by the Boltzmann-Gibbs 
statistics. The Boltzmann-Gibbs entropy is not proportional to the number of particles Ar in the reservoir and 

is not extensive. The extensive quantity is the generalized entropy Sq ' for some value of g 7 ^ 1: Sq ~ Ar. Here 

1 

we consider only the situation when <7 < I. The number of microstates in such a reservoir is Wr = 65 ®’ ~ ■ 

Therefore, the Boltzmann-Gibbs entropy = k^ In Wr depends on the number of particles in the reservoir as 
5'r In Ar. This expression for the Boltzmann-Gibbs entropy is similar to the entropy of the system consisting 

of d = quasi-particles, whereas the number of particles Ar plays the role of the volume. A simple model of such 
a system has been presented in Ref. [ 3 ^ : the model consist of a spin chain containing Ar spins; spins next to each 
other have almost always the same direction, except there are d cases when the next spin has an opposite direction. 
In the Boltzmann-Gibbs statistic large reservoir has large heat capacity Cr. Similarly, here we require that the < 7 - heat 
capacity of the reservoir, defined by Eq. (1551) . should be large. 

The approach presented in this Section is similar to the approach in Ref. [3l|. However, in Ref. [3l| the reservoir is 
incorrectly interpreted as a heath bath and having large heat capacity. As we have seen in the previous Section, such 
a reservoir leads to the exponential Boltzmann factor and approximation of the expansion of the number of states 
as a q-exponential is not justified. Interaction of the system S weakly coupled to a finite reservoir having a finite 
energy has been considered in Ref. [d^. Under the assumption that the number of microstates of the reservoir having 
energy less than Ar grows as a power-law of Ar, the q-exponential distribution of the energy of the system has been 
obtained. However, in Ref. [d^ the extensivity of the generalized entropy has been not used and the parameter q 
tends to I when the number of particles of the reservoir increases. 

As in the previous Section, the system under consideration S is interacting with the reservoir via short-range 
interactions, thus the total energy of the composite system in the macroscopic limit is Atot = E + Ar and the 
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total number of microstates in the composite system when the system S has energy E is W{E)W’B,{Etot — E). 
Introducing the generalized entropy of the system Sq{E) = In^ W{E) and the generalized entropy of the reservoir 

Sq^\Eji) = fce In, WR(illR) we can write 


WtotiEtot) 


E 


T^S,{E) 

Cq Cq 


E 


(32) 


When g < 1 this sum cannot be approximated by the largest term. Approximation of a sum of large q-exponentials 
is investigated in Appendix We assume that the postulate of equal probabilities of microstates in the equilibrium 
is valid also in the non-extensive statistical mechanics. When the postulate of equal probabilities of microstates is 
assumed, the statistics of the system S is determined by the reservoir according to Eq. m- Therefore, even an 
ordinary system interacting with the reservoir having large q-heat capacity is more conveniently described by the 
q-entropy. 

The largest term in the sum (I32p corresponds to the most probable state of the composite system and is found from 
the equation 


AS,{U) sfc4"’(E.o.-Cf) 

Here U is the most-probable energy of the system S. Thus in order to satisfy the zeroth law of thermodynamics we 
need to define the temperature T as 

1 S,S,iU) 

T 1 + _ (/) ■ 

This definition of the temperature is the same as Eq. m- 

If one introduces the entropy of the combined system as = A:b In^ lE(Etot) then as a consequence of 

the impossibility to approximate the sum (I32|) by the largest term the entropy of the combined system is not a 
simple combination of the entropies of the system S and the reservoir R: Sq^°^\Etot) ^ Sq{U) + Sq^\Etot — U) + 

^^Sq{U)Sq^\Etot — U). Due to this the conclusions of Ref. [dlj that the zeroth law of thermodynamics holds only 
if the energy is also nonadditive does not apply for the situation considered in this paper. On the other hand, if 
one assumes that the interaction between the system S and R are long range and the energy is not additive then the 
pseudo-additivity of entropies can be valid [s^ . 

The probability of a microstate of the system S is given by Eq. ®. Similarly as in the previous Section, assuming 
that the second derivative of g-entropy of the reservoir is very small, Sq^\Etot) ~ 0, we can approximate 


WniEtot -E) = e; 




Sf) (Btot -tf) -^(E-U) ^ Sf) (Btot -U) 


(34) 


Note, that now the Boltzmann-Gibbs entropy of the reservoir does not have a small second derivative. The condition 
at'i Sq^\Etot) ~ 0 means that the g-heat capacity of the reservoir, defined in Eq. (l38|) . is very large. The ordinary 
heat capacity of such a reservoir can be determined as follows: if we increase the energy of the reservoir by a small 
amount AE, the increase of the temperature T, according to Eq. (1551) . is 


AT = 


1-g 

kB 


AE-T- 


q(R) / \ 

[^*9 \ Eb ,) 

d c(R) / ^ \ 
OEr^I (Ar) 


AE. 


This means that the heat capacity of the reservoir Cr = is 


Cr 


1 


1-9 I T 1 ’ 

fcB 


where Tq^'^ 


is the auxiliary g-temperature of the reservoir defined via the equation 


1 

ry-,(R) 


d 




(35) 


(36) 


( 37 ) 
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and 


^(R) ^ 

9 




( 38 ) 


is the g-heat capacity of the reservoir, defined similarly to the physical heat capacity, Eq. (1^ . If the second derivative 

o2 f'R') U 

of the generalized entropy is small,'(E r) Ri 0, then the heat capacity is C = If we increase the energy of 
the reservoir with very large g-heat capacity by AE, the new temperature of the reservoir becomes 


T' = T+ ^—^AE. 
Kb 


(39) 


The expression for the heat capacity C = is the same as the heat capacity of a gas consisting of d = 
quasi-particles. 

Using the property (IB2I) of the g-exponential function we obtain that the probability of the microstate of the system 
S is proportional to the factor 


where T is the temperature according to Eq. (1331) . However, the temperature T depends also on the properties of the 
system, not only on the reservoir. It is more convenient to introduce the temperature that the reservoir not interacting 
with the system could have: 


1 ^ 

Taking into account that 0 we obtain 


T{U) 


1 I 1—9c(R)/'z 7' ^ 1—9rr d c 

^ + SETT*' 


fce 9Et, 


dR) 

^9 


d c(I^) 

OEtot ‘^<1 


(E„) 




k-B 


(41) 


(42) 


This equation shows that the interaction with the system lowers the temperature of the reservoir. On the other hand, 
the g-temperature of the reservoir, defined by Eq. dSH) does not change. 

Using Eq. (1421) we get that the probability of the microstate of the system S is proportional to 


. (43) 

An expression similar to Eq. (H51) has been obtained in Ref. [3l|. However, in Ref. [3l| the temperature that enters 
P{E) has been interpreted as a physical temperature T, because the reservoir has been assumed to be a thermostat. 
The correct observation that the energy of the reservoir interacting with the system should decrease has been presented 
in Ref. [H]. 

One common objection to Eq. (H51) is that this expression is not invariant to the change of zero of energies [I|. 
However, this reflects the physical situation of the system interacting with the reservoir having very large g-heat 
capacity and, consequently, small physical heat capacity. The zero of the energy of the system S is fixed by the 
requirement that the energy of the reservoir should be Utot when U = 0. If we shift the energy zero by AE, this 
is equivalent to the decrease of the energy of the reservoir by AE. This decrease of the energy of the reservoir 
decreases the temperature. The probability of the microstate should remain the same, thus the new factor should be 
proportional to the old: 


P'{E) = expg 


fcBT'(O) 


E 


P{E + AE) = exp - 


A:bT(0) 


{E + AE)] . 


(44) 


It follows that 


1-g 

/CB 


AE. 


T'(0) = T(0) 


(45) 
















This equation is consistent with Eq. (1391) . Similar argument has been presented in Ref. (40j | by considering a system 
S interacting with a finite reservoir. 

Interesting feature of Eq. (l43ll is the presence of the cut-off energy: it follows from the definition of the (/-exponential 
function that P{E) becomes zero when E ^ Emax where 


i^max = -^T(O) . 
l-q 


(46) 


This property of P{E) ensures that the physical temperature T is always positive. Discussion of possible cut-off 
prescriptions associated with Tsallis’ distributions is presented in Ref. (dl j. 

Using the factor (l43l) we can write the normalized probability of the microstate as 


P(K) = 


(47) 


where 




"fcBr’Co)'®'* 


(48) 


is the generalized partition function. The distribution of the energy of the system E is obtained multiplying the 


probability p{^i) by the number W {E) = el 


of microstates having energy E^ = E: 


P{E) = —eq 

^q 




- y eq, 


where 


T(£;) = T(o) - ■ 

Kb 

The probability p{E) should be normalized, thus the partition function can be also written as 




1 T (E) q (p\_ 
fcg T(0) 


kBT{0) ^ 


E 


The energy U corresponding to the largest term in the sum (151|) is determined by the equation 

mfS.iU) 1 


1 + T{u) • 

According to Eq. the temperature T{U) coincides with the physical temperature. 


(49) 


(50) 


(51) 


(52) 


IV. GENERALIZED FREE ENERGY 

The probability proportional to the factor (H51) admits several different possibilities to generalize the free energy. 
First of all, there are three possibilities corresponding to three temperatures: initial temperature of the reservoir 
T(0), auxiliary (/-temperature Tq and the physical temperature T{U). From those three choices only the temperature 
T(0) depends only on the reservoir and does not depend on the properties of the system. On the other hand, the 
temperature T(U) has a direct thermodynamical interpretation. In addition, the average energy of the system is 
connected to the generalized entropy with the parameter 2 — q. The derivative of this generalized entropy with respect 
to average energy yields another auxiliary temperature T 2 -q and the corresponding generalized free energy. 


A. Initial temperature of the reservoir and unnormalized (/-averages 

Let us consider first the generalized free energy Eq corresponding to the temperature T(0). This choice is closely 
related to the approximation of the sum of large (/-exponentials and to unnormalized (/-averages. When q < 1 the sum 
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in Eq. (ED cannot be approximated by the largest term even in the macroscopic limit. The approximate expression 
for the sum of large g-exponentials is obtained in Appendix According to the results of Appendix and Eq. (ISTl) 
the (/-logarithm of Zq can be approximated as 


For any function of the energy /(A) the following equality holds: 

1 / 1 a —1 


E E,. 


^Sq{E^) / ■^Sq{Efj,) 




f{E^) 


l+^-^Sq{E,) 




Therefore, we can approximate the (/-logarithm of Zq as 

^SqjE,) _ 

^l+l^Sq{E^) WO) 


1 Pipr . 


^ \ “U / 

This equation suggest to introduce the unnormalized (/-average energy of the system 

Uq = Y.EWY = Y,E 


l + K^Sq{E)]pW 


fce 


(54) 


(55) 


(56) 


This unnormalized (/-average of the energy can 
equation 

Uq = fcBr(0)2 


be determined from the generalized partition function Zq using the 

^ ’ " (57) 


dT{0) 


ln„ Z, 


l(j ■ 


In analogy to Eq. ED we introduce the generalized free energy 

Eq = -kBT{0)\nqZq. 

The equation 

Fq=Uq- T{0)Sq 

defines the entropy Sq which is related to the unnormalized (/-average of 
obtain 

dSq ^ 1 

dUq T(0) ■ 


(58) 


(59) 


the entropy Sq. Using Eqs. ED^dSD we 

(60) 


Entropy Sq can be calculated using the probabilities p{p) according to Eq. dD. Indeed, we have 

Sq = 7^.iUq- Eq) = ^ T^^E S ■ 


(61) 

(62) 


Expressing the energy from the probability p{p) we get 

Efj. = -kBT{0)\a.q[p{p)Zq\ = -kBT{0) (lngp(/x) + p{p)^~'^ Iwq Zq) . 

Inserting this expression for the energy E^ into Eq. ED and taking into account that J2fj,Pip) = 1 we obtain 

- 1 


Sq — ■ 


“ 1-q 

This expression is consistent with the approximation (1551) . According to the approximation (I53p the entropy Sq is 


.y:s,(E)p(E)« = y:—|t|— 




In the macroscopic limit the entropy Sq{E) is large and we can approximate 

E„p(c)« -1 


1 — (/ 


■kB^^^ 


1-9 


(64) 


(65) 


This expression is the same as ED- 
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B. 5-temperature and normalized 5-averages 


It can be more convenient to deal with normalized 5 -averages. The normalized 5 -average of the energy is 

^ ■ 

Using Eq. (El the normalized 5 -average of the energy can be written as 

^ _ EEE{l+^-^S,iE))piE)<^ 

EE{^+^-^S,iE))piE)<^ 


( 66 ) 


(67) 


The sums of the form f{E)p{EY can be approximated as f{U) '^eP(E)‘^ + f'{U) J2EiE — U)p{EY . The sum 
J2e(E — U)p{EY is small, since close to the maximum U the probability p{E) is an even function oi E — U. We 
obtain that in the macroscopic limit the normalized 5 -average of energy Uq should be close to the most probable 
energy U. Since the sum of large 5 -exponentials cannot be approximated by the largest term when 5 < 1, it is not 
possible to determine the most probable energy U or the average energy Ef^p{p) knowing only the generalized 
partition function Zq. However, it is possible to calculate Uq, which is close to U. On the other hand, the entropy Sq 
cannot be approximated by Sq(U). 

Using Eq. (1631) the normalized 5 -average Uq can be related to the unnormalized Uq via the equation 


Uq= - 

^ + T^Sq 


( 68 ) 


From Eq. (16811 follows that the introduction of the normalized 5 -average energy Uq allows to factorize the generalized 
partition function Zq-. 


1 q _ 1 TT 1 q _ 1 TT 

y _ fcBT(O) fcBT(O) ^<2 

Zjq — tq — fiq t-q 


(69) 


The entropy of the combined system using the expansion (1341) can be written as 


^dot) ^ ^\f^^Zq 


(70) 


Using Eq. (l69)l and assuming that U fnUq we get that the entropy Sq^°^'^ can be expressed as the usual pseudo-additive 
combination of entropies from the non-extensive statistical mechanics: 


^(tot) ^ S^^^Etot - Uq) + Sq + HUtot - Uq)Sq . 


(71) 


However, in this equation the entropy Sq is not directly connected to the number of microstates of the system S. 
Let us introduce an auxiliary 5 -temperature Tq of the system S via the equation 


Using Eqs. EQl) and (| 68 )l we get 


1 dSq 
Tq dUq ■ 


(72) 


TiUq) = Tq(^l + ^Sq'^ , ( 73 ) 

where T{Uq) = T(0) — ^^Uq is the temperature of the reservoir corresponding to the energy of the system equal 
to Uq. Since Sq > 0, the 5 -temperature is always smaller than the physical temperature T(Uq). Note, that only 
physical temperatures of the system and the reservoir are equal. The 5 -temperature of the system Tq is not equal to 
the 5 -temperature of the reservoir Tq^'^. 

We introduce the 5 -analog of the free energy corresponding to the temperature Tq-. 

Fq=Uq- TqSq . 


(74) 
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Then, using Eqs. (EH) and (Ell) we get 


’ dT^ ■ 


We define the g-heat capacity of the system as 


dU, dS, d^F, 

dTq ‘‘dTq dry 


( 75 ) 


(76) 


The physical heat capacity C can be determined as the derivative of Uq with respect to the physical temperature 
T{Uq): 


C= 


(77) 


dT{Uq) • 

Using Eqs. El and El we get the equation that relates the physical heat capacity with the auxiliary g-heat capacity: 

1 


C = 


T{U^) 1 , ■ 

T, C, fcB 


(78) 


Since Tq < T{Uq), from Eq. ((75)) follows that the physical heat capacity C is always smaller than the g-heat capacity 




Generalized partition function Zq related to the generalized free energy Fq is 


Z = p ^ “ 

Z/q — Pq 


Note that Zq^ Zq. Using Eqs. El, (El and El we get the expression for the energy Uq: 

Uq = k^Tq In^ Zq . 


(79) 


(80) 


The generalized partition function Zq cannot be directly expressed as a sum. However, Zq can be connected to a sum 
of g-exponentials as follows: we write the probability of the microstate in the form 


PW = 


(81) 


where 




“ k^T{U,) ^l) 


(82) 


is related to Zq via the equation 


Using Eq. (l69ll we get 


Zq = ZqCq 


1 TT 

tTtTTX ^ a 


(83) 




(84) 


Therefore, 


lug Zq = In, Zq - jr^Ug . 


( 85 ) 
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C. Physical temperature and Renyi entropy 


The third possibility is to in introduce the free energy corresponding to the physical temperature T(Uq). In order 
to do this let us consider another entropy, given by the equation 


Sq = k-Q In eq 


— \n(l 
-<1 \ 



( 86 ) 


The entropy Sq is more directly connected to the physical temperature T{Uq). Indeed, using Eqs. (1721) and (|73)) we 
get that the derivative of the entropy Sq gives the physical temperature: 


From Eq. (ESI) it follows that 


dSq ^ 1 

dUq T{Uq) ■ 




(87) 


( 88 ) 


Thus the entropy Sq is the Renyi entropy [H, Q . 

Using the Renyi entropy we introduce the free energy corresponding to the physical temperature T(Uq): 

Fq = Uq- T{Uq)Sq . (89) 

Also in this case we retain the Legendre transformation structure. For example, using Eq. EH) we get 


dPq 

dT{Uq) 


(90) 


D. Average energy of the system 
It is impossible to exactly determine the average energy of the system 

U = E 


(91) 


knowing only the sum of g-exponents Zq. However, the knowledge of another sum 


E 




2-9 


allows us to do so. Indeed, using the property of the ^-exponential function (IB4|) and the expression for the probability 
p{^j.) Eq. (ITtI) we get 


- _ fcB'r(o)^ d ^ 

“ (2 - q)Zq 5T(0) 4- 


'ksT{0)^l^ 


2-9 


(92) 


Instead of this sum we can use the generalized entropy o with the parameter q' = 2 — q: 


S2-q — k-Q- 


1-9 


Using the generalized entropy S' 2-9 the expression for the average energy U becomes 

- _fcBlW 2-9 f. l-gp 
(2-q)Z,9r(0)^« fcB 


(93) 


( 94 ) 
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We can obtain another expression for the average energy U by inserting from Eq. (l62l) into Eq. (IM]) : 

U = T{0)Zl-‘>S2-g - kBT{0) In, Z, . 

Combining Eqs. (|M|) and dSH) we obtain 

du r(o)zi-9 dS 2 -g 

dT(0) 2-q dT{0) 

Therefore, yet another auxiliary temperature T 2 _,, introduced by the equation 


(95) 


(96) 


1 _ dS2-g 

T2-g dU 


is equal to 


Z^-g 


(97) 


(98) 


This relation between temperatures is exactly the same as obtained by maximizing the entropy (1931) with the constraint 
m 0- The temperature T 2 _, , similarly as the temperature T,, depends not only on the reservoir but also on the 
properties of the system. 

The generalized free energy corresponding to the average internal energy of the system U and the temperature T 2 _, 
is 


F2-q = U- T2-gS2-g • (99) 

This expression for the generalized free energy E 2 _, is similar to the expression 

F, = C7 - (2 - q)T2-gS2-g (100) 

for the generalized free energy F, that follows from Eq. (IM|) . We see that in general F, ^ F 2 _, . 


V. DISCUSSION 

In summary, we have demonstrated that a small system interacting with a large reservoir having large g-heat capacity 
can be described by the non-extensive statistical mechanics. From the point of view of the ordinary statistics such a 
reservoir is similar to a gas oi d = l/{l — q) quasi-particles. The probability of the microstate of the system interacting 
with the reservoir via short-range forces is given by the g-exponential function (14311 . instead of the exponential 
Boltzmann factor (nnn- Large g-heat capacity of the reservoir leads to a small physical heat capacity, therefore the 
temperature in the equilibrium T depends both on the properties of the reservoir and the properties of the system. 

In order to avoid this inconvenience one can consider the temperature T(0) of the reservoir that is not interacting 

f R,^ 

with the system or introduce an auxiliary g-temperature T, ' that remains constant due to large g-heat capacity of 
the reservoir. Small heat capacity of the reservoir does not allow to consider it as a thermostat, thus the description 
using the standard canonical ensemble of the statistical mechanics is not applicable. The treatment of the canonical 
ensemble presented in this paper allows us to obtain relations between the physical temperature T and the auxiliary 
g-temperature F, (17^ as well as between the g-heat capacity C, and the physical heat capacity C (1781) . 

Sums of large exponentials often appear in the Boltzmann-Gibbs statistical mechanics. Such sums can be approx¬ 
imated by keeping only the largest term. Similarly, in the non-extensive statistical mechanics appear sums of large 
g-exponentials. However, for such sums taking only the largest term is a very poor approximation. This is because 
the g-exponential function with g < 1 does not decrease as fast as the exponential function. As a consequence, the 
deviations from the most probable state in the non-extensive statistical mechanics are much larger than the deviations 
in the standard statistical mechanics. As it is shown in Appendix [B] sums of large g-exponentials are well approxi¬ 
mated using g-averages. This fact is one the reasons why g-averages play such an important role in the non-extensive 
statistics. 

In this paper we considered the reservoir for which the generalized entropy with g < 1 is extensive. The case of 
g > 1 is more complicated, because one cannot take the macroscopic limit N —>■ oo. The investigation of the small 
system interacting with the reservoir characterized by g > 1 remains a task for the future. 
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Appendix A: Sum of large g-exponentials 


One can easily show that the sum of large exponentials 

IV 

Z = J2 

can be approximated by the largest term. Indeed, if (/)max is the maximum of then 

gAf^max ^ ^ ^ p{/gAf0max 


and 


0 € 


InZ 




InW 

N 


(Al) 


(A2) 

(A3) 


If W grows slower than exponentially with increasing N, then in the limit of large N the ratio InW/N vanishes and 
we have 


lim 


InZ 


N^qo ]V 

Now let us consider the sum of large g-exponentials 

w 




>N(p{i) 


i=l 


(A4) 


(A5) 


where q < 1. In contrast to the sum of large exponentials, approximation of the sum of large g-exponentials with the 
largest term is a poor one. We can construct a better approximation as follows: let us introduce the weights 

N4>(i) 

Pi^) = ^ (A6) 

and the unnormalized g-average 

(</>), (A7) 

i-1 


By noticing that 

N(j){i) = \nq[p{i)Zq] = Inq p{i) + In, Z, 

one can write the difference In, Z, — N{(j))q as 


ln,Z,-iV(</)), = ^i^lh!)LJ 


(AS) 


(A9) 


The sum have the largest possible value when all weights p{i) are equal. In such a case p{i) = XjW and 

Y,,p{iY = Thus 




N 


N 


If W grows with increasing N as and 


p < 


l-q 


then in the limit of large N the ratio In, W/N vanishes and we have 

lim = {4>)q ■ 

Af->oo N 


(AlO) 

(All) 

(A12) 


This gives the required approximation of the sum of g-exponentials. Note, that In, W/N as the upper limit of the 
difference In, Z, — N{(f>)q is the worst case, when all terms in the sum are equal. For sufficiently fast decreasing terms 
the sum be bounded even for larger W. 
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Appendix B: Some properties of g-exponential function 


In this paper we have used the following properties of g-exponential and g-logarithm: multiplication of two q- 
exponentials 


^[l + {l-q)v]x+y _ ^x+[l + (l-q)x]y 



(Bl) 

(B2) 


p^+V 


g-logarithm of a product 

In, xy =[l + {I - q) In, y] In, a: + In, j/ = In^ x + [1 + (1 - q) In, x] In, y , 
the derivatives of g-exponential and g-logarithm: 


(B3) 



(B4) 


(B5) 


The ^nations (IB1I) ~ (IB5I) can be easily derived using the definitions ([5]) and are presented in the Appendix A of 
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